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ABSTRACT 

The elasticity problem for a long hollow circular cylinder contain- 
ing an axi symmetric circunferential crack subjected to general nonaxi sym- 
metric external loads is considered. The problem is formulated in terms 
of a system of singular integral equations with the Fourier coefficients 
of the derivative of the crack surface displacement as density functions. 
The stress intensity factors and the crack opening displacement are cal- 
culated for a cylinder under uniform tension, bending by end couples, and 
self-equilibrating residual stresses. 


1. Introduction 

The elasticity problem for a cylindrical structure such as a pressure 
vessel or. a pipe which contains a part- through surface crack appears to be 
analytically intractable. Such problems are generally treated either 
numerically by using the technique of the finite elements [1,2] or the 
boundary integral equations [3], or, in relatively thin-walled cylinders, 
approximately by using the line spring model in conjunction with the shell 
theory [4,5]. The limiting cases of some of these problems can also be 
solved analytically which provide very useful results, for example, regard- 
ing the bounds of stress intensity factors. The plane strain problem of 
a hollow cylinder containing a crack in a radial plane which constitutes 
the limiting case of a cylinder with an axial part-through crack is one 
such problem. Another problem is that of a long thick-walled cylinder 
which contains an axisymmetric radial crack. In this paper the latter 

(*) This work was supported by NSF under the Grant CME-78-09737, by 
NASA-Langley under the Grant NGR 39-007-011 and by DOT under the 
contract D0T-RC-82007. 
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problem is considered for arbitrary nonaxi symmetric loading conditions. 
The related problem for an elastic solid cylinder containing a penny- 
shaped crack is considered in [6-8], and that for a hollow cylinder under 
axisymmetric loading conditions in [9]. 


2. Formulation of the Problem 

The problem under consideration is described in Fig. 1. The external 
loads may be decomposed in such a way that the problem may be expressed 
as the superposition of a problem which is symmetric and that which is 
anti -symmetric with respect to the z=0 plane. In this paper, only the 
symmetric problem is treated. Clearly, it is sufficient to consider one 
half (e.g., z>0) of the medium only. Furthermore, for the given external 
loads the quasi static problem for the cylinder without the crack is 
assumed to have been solved. Thus, the information regarding the stress 
intensity factors may be obtained by considering the perturbation prob- 
lem in which the crack surface tractions are the only external loads. 
Following are then the boundary and continuity conditions of the problem: 


tf rr (a,0,z) = T rz (a,e,z) = T re (a,e,z) = 0, 0<e<27r, 0<z<» , (1) 
a rr (b,e,z) = T rz (b,e,z) = T r0 (b,0,z) = 0, O<_0<2t it, 0<6<2Tr, (2) 
T zr (r, 0 ,O) = T z0 (r, 0 ,O) = 0, a<r<b, O<0<2 tt , (3) 
cr zz (r,0,O) = p(r,e) , c<r<d , 0<e<2ir , (4a) 
u z (r,e,0) = 0, a <r<c, d<r<b, 0<e<2ir , (4b) 


where the dimensions of the cylinder and the crack are given in Fig. 1 
and p(r,e) is a known function. Also, the stresses must vanish as z-*». 

Referring to [10], [11], and [8] the solution of the problem may 
be expressed in terms of a system of five harmonic functions as follows: 
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The harmonic function 4 > is associated with the formulation of the 
problem for an infinite elastic space for which z=0 is. a plane of symmetry 
[10]. The functions B 0 , Bp B 2 and $ are equivalent to Papkovich-Neuber 
potentials in cylindrical coordinates [11,8]. With an important application 
of the bending of the cylinder in mind, if we restrict the considerations 
to external loads which are symmetric in e, that is, if p(r,e) = p(r,-e), 
then the functions p and <j> may be expressed as 

(18) 

(19) 

where the functions A n , (n=0,l,...) are unknown and 

ir ir 

o 0 (r) = p(r,e)de, a n (r) = \ J p(r,e)cos node, (n=l ,2, — ) . 

o 0 

( 20 ) 

For the symmetric loading under consideration from (5)-(17) it may 
be seen that B Q and B-j are even and B 2 and ip are odd in e. Taking also 
into account the symmetry with respect to the z=0 plane, these harmonic 
functions may be written in the following form: 


p(r,e) = z a (r) cos ne , 
o n 


00 f —rt7 

4 >(r,e, z) = z cos ne A n (a)J n (ctr)ae da, 
0 J 
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00 

B 1 = „f 0 cosCCn+l)e3 f | [C ln (s)I n+1 (sr) + C 2n ( S )K n+] (sr)]cos(sz)ds, 

0 

( 21 ) 

00 09 

B 2 = r sin[(n+l)e] | | [C ln (s)I n+1 (sr)+C 2n (s)K n+1 (sr)]cos(sz)ds , 

( 22 ) 

00 

B 0 = J Q cos(ne) || [C 3r (s)I n (sr)K: 4n (s)K n (sr)]cos(sz)ds , (23) 

00 

\p = z sin(ne) ! [ C C 5n<^ n( sr ) +c 6rf s )y sT ^ cos ( sz ) ds • (24) 

n=0 J o 


where I n and K n are modified Bessel functions and the functions C in , 
(i=1,...,6; n=0,l,...) are as yet unknown. 

From the formulation of the problem as stated, it may be verified that 
the conditions (3) are identically satisfied. The seven sets of unknown 
functions A n , C-j n ,..., Cg n , n=0,l,... are then obtained from the seven 
conditions (1), (2), and (4) by observing that these conditions are in the 
form of sine or cosine series, and hence, by writing the coefficients of 
sin ne or cos ne equal to zero for each n. Thus, for each n the homogeneous 
conditions (1) and (2) are used to eliminate six of the unknowns and the 
mixed boundary condition (4) is used to determine the remaining unknown 
function. The integral equation of the problem may be obtained directly 
by expressing 

u (r,e,0) = z <j> n (r)cos ne, a£r<b, 0£6<2ir, (25) 

where from (4b) it follows that 

<f> n (r) = a<r<c, d<r<b, . (n=0,l, ..'.,) (26) 
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From (7), (19) and (21)-(26) it follows that 
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By using the conditions (1) and (2), the unknown functions C-. 
C 6n may be obtalne d in terms of <}> n in the following form: 





d 

m ij( s )j 
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(28) 


where the functions m^- and Gj, (i,j=l . ,6) are given in Appendix A. 
Finally, from (4a) and (18), it follows that 


z-Io -0 + az)A n (a)J n (ctr) a 3 e aZ da 
0 

oo 

+ \ | S E C in (s)N i (r, s)ds = a n (r), (n=0,l ,.. .) , (29) 

o 11 

where 

N-j = 2vl n (sr) + sr I n+1 (sr), N 2 = -2v ^(sr) + sr K n+1 (sr), 

N 3 = s I n (sr), = s ^(s r) . (30) 

Substituting now from (27) into (29), through an asymptotic analysis 
separating the singular part of the kernel, and by integrating by parts 
we obtain 
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where the kernels L-j and Lg are given by 
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The function ^(z) which appears in (29) is the Lommel's function 
[12,13], and C(x) and S(x) are the Fresnel integrals defined by [14] 


X 

C(x) = -L f -SSU. dt , S(x) • -L f dt . (34) 

/2tr * q /t /2ir ■'q /tf 

The functions v-(s,t) appearing in (33) are given in Appendix B. A conven 

J 

ient numerical technique for the evaluation of Lommel’s functions of real 
and imaginary arguments is described in Appendices C and B. 
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The infinite integrals which appear in (32) converge rather slowly for 
nearly all values of r and t. The reason for this is that, considered 
term by term, all integrals in (32) are divergent. However, it can be 
shown that the asymptotic behavior of the terms for large values of a is 
identical but opposite in sign to the trigonometric terms in the integrand 
and the sun leads to a convergent integral. The real ted asymptotic analy- 
sis and the numerical evaluation of these integrals are described in Appen- 
dix D. 

3. The Stress Intensity Factors 

The Mode I stress intensity factors along the crack borders r=c and 
r=d are defined by 

k(c) = lim ^(c-r) a (r,e,0) , (35) 

r-K: 

k(d) = lim /2(r-d) 0__(r,e,O) , (36) 

r-Ki zz ■ 

and, in terms of the crack surface displacement, may be expressed as 

He) = lim /2(r-c) ^ u (r,e,0) 
r->c 

= yL lim /2(r-c) z ^ (r) cos ne (37) 

r->c o 

k(d) = - rh lim A- u (r,o,0) . 

1 v r-»d 3r z 

= - Tallin /2(d-r) z <L(r)cos ne . (38) 

r-Ki 0 n 

The formulation given in the previous section is valid provided the 
external load p(r,e) is such that the crack surface displacement u z (r, 8,0) 
is positive everywhere in (c<r<d, JO^u&r), — -In— this case, the stress inten- 
sity factors k(c) and k(d) would be positive for all values of e. 
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4. The Embedded Crack 

Referring to Fig. 1, if c>a and d<b the crack is embedded and the crack 
surface displacement is zero at r=c and r=d. From (25) it then follows 
that 4 > n (c ) = 0 = <j> n (d), the integrated terms on the left-hand side of the 
integral equation (31) vanish and the integral equations must be solved 
under the following conditions (see (26): 
d 

| <i>p(t)dt = 0 , n = 0,1,... (39) 

c 

The solution of the singular integral equation (31) is of the form 

♦„(t) = f n (t)[(t-c)(d-t)] -i , c < t < d , (40) 

where f n (t) is a bounded function. After normalizing the interval (c,d) by 
defining 

t = , r = , f n (t) = F n (0, (41) 

The functions F n (-r), n=0,l,... may be obtained from (31) by using a Gauss- 
Chebyshev type quadrature formula [15]. 

5. The Edge Cracks 

If the crack is a surface or an edge crack, then the crack surface 
displacement will be zero only at one tip and consequently the conditions 
(39) will no longer be valid. In this case (the Fourier coefficients of) 
the crack opening displacement on the cylinder surface may be expressed as 

d 

a = c < d < b: <|> n (d) = 0 , <j> n (c) = -J <^(t)dt , (42) 

c 

f d 

a < c < d = b: 4> n (c) = 0 , <}> n (d} = j <j>^(t)dt , (43) 

c 
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for the inner and the outer edge cracks, respectively. For these crack 
configurations, even though the integral equation (31) is still valid, the 
integrated terms on the left-hand side are not zero and by using (42) 
and (43), (31) may be expressed as 


'<![(£> (f) ♦ LjCt.r) - L l( s,r) 


+ L 2 (t,r)-L 2 (s,r)}<f>'(t)dt = ^-o n (r), c<r<d, n-0,1,... (44) 


where s = c = a<d<b for the internal edge crack and s = d = b>c>a for 
the external edge crack. . 

For the edge cracks even though the kernel L,(t,r) contains addi- 
tional singular terms (in the form of generalized Cauchy singularities) 
which become unbounded as t and r go to the end point s simultaneously, as 
seen from (44) for (t,r) s the entire kernel in (44) vanishes. As 
shown in [16], the consequence of this is that the derivative of the 
crack surface displacement at r=s (i.e., on the cylinder surface) is 
bounded. This point is rather important in the numerical solution of 
the integral equation. For the embedded crack, in the Gaussian integration, 
referring to (41) the unknowns and the collocation points are F n ( t.- ) and 
p., respectively where 

J 

Tj = cos (it jjjpp) , i = 1 m , (45) 

Pj . = cos (it |j^-) , j = l,...,m-l . (46) 


The additional equation to solve for m unknowns F n (r i ), (i=l,...,m) is 
provided by (39). However, in the edge crack problem (39) is not valid 
and consequently there are only m-1 algebraic equations. In this case 
since the column in the coefficient matrix obtained from (44) which cor- 
responds to t=s, i.e., t=1 or t=-1, is identically zero, the system of 
algebraic equations involve m-1 unknowns only and hence gives a unique 
solution. 
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After obtaining the stress intensity factor at the crack tip may 
again be obtained by using (37) or (38). In the edge crack problem 
another quantity of some interest is the crack opening displacement at the 
cylinder surface, which, referring to (25), (42), and (43), may be 
obtained from 


d 


s = u (a,e,+0) - u (a, 0 ,- 0 ) = -2s cos ne 

z z n J 

1 

<J>p(t)dt , 

(46) 

a 


b 

0 


s = u (b,e,+0) - u (b,e,-0) = 2s cos no 

■ z 0 J 

4»„(t)dt . 

(47) 


6. Results 

The numerical results given in this section are obtained for three 
different loading conditions, namely the uniform axial stress, the bending 
of the cylinder by end couples, and the self-equilibrating residual stress. 
For these three loading conditions the crack surface tractions in the per- 
turbation problem are respectively given by (see (4a) and (18)) 

p(r,e) = a Q (r) = . (48) 

p(r,e) = O^rjcose = -^(j^cose , <j-| = (49) 

P(r,e) = « 0 (r) = -„ s C ^g l ^rl - 1] , (50) 

where P*, is the axial force, M is the bending moment, and a s is the mag- 

nitude of the compressive stress on the surfaces of the cylinder. In 
the residual stress problem the axial stress a z (r,e,0) is assumed to be 
independent of e, parabolic in r, compressive on and near the surfaces, 
tensile in the interior region of the cylinder wall, and statically self- 
equilibrating. Unless otherwise stated, in all examples the Poisson's 
ratio is assumed to be 0.3. 
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For a ring-shaped crack in an infinite solid (i.e,, for a=0, b=®) 
the calculated stress intensity factors are shown in Figure 2. Here a Q 
is the uniform crack surface pressure, for n=l the crack surface traction 
is defined by 

p(r, 0 ) = CT-j(r)cose = -a-j (^j-)cose , (51) 

and ^ 

k = k(s,e) = k'(s) n cos ne , (n=0,l) (52) 

is the stress intensity factor along the crack edges s=c and s=d. Note 
that as c/d approaches 1 and 0, respectively, the plane strain and the 
penny-shaped crack results are recovered. 

For hollow cylinders under tension or bending the results are given 
in Tables 1-17. Tables 1-3 show the stress intensity factors for a sym- 
metrically located imbedded crack (i.e. , for b-d = c-a. Figure 1). Here 
the normalizing stresses a Q and 0 -| are defined by (48) and (49). As 
(a/b) 1 the curvature effect disappears and the stress intensity factors 

approach those obtained from the plane strain solution of a strip contain- 
ing a center crack. 

Tables 4-15 show the stress intensity factors and the crack opening 
displacements for an internal or an external surface crack. Here Z 
is the total crack depth ( Z = d-a for c = a and Z = b-c for d=b), h = 
b-a is the wall thickness (Figure 1), and the crack opening displacement 
6 is defined by (46) and (47). For a very small crack depth (i.e., for 
Z/h - 0.01) the effect of a/b on k is shown in Figure 3. In this case 
the stress intensity factors for a half plane, for a strip and for a 
solid cylinder are practically the same (i.e., k /a Q /£ = 1,121). It is 
seen that for the internal edge crack as a/b approaches 0 and 1 k/a Q vT 
approaches respectively 2/tt and 1.121 which are the values for a penny- 
shaped crack and for a strip. For the external edge crack, on the other 
hand the influence of a/b on k is hardly noticeable. 

The limiting values of k and 5 for (a/b) 1 shown in the tables 

are obtained from the plane strain solution of a strip containing an 
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edge crack. The stress intensity factor for an internal edge crack seems 
to be a monotonously increasing function of a/b for all crack depths. On 
the other hand, for an external edge crack k always seems to go through 
a minimum as a/b goes from 0 to 1. 

For one crack-cylinder geometry, namely for a=c<d<b, (a/b) =0.5 and 
(£/ h) = 0.3 the effect of Poisson's ratio on k and s is shown in Tables 
16 and 17. It may be seen that both k and 6 increase monotonously but 
very slightly with increasing Poisson's ratio. 

For the residual stress problem the crack geometry and the stress 
profile are shown in Figure 4. It is clear that if the crack is suffi- 
ciently shallow so that it lies in the compressive zones near the surfaces, 
then the crack surfaces would remain closed and k would be zero. In the 
example under consideration by letting a Q (r .) = 0 (i = 1,2) in (50) we 
find 


r.-a r^-a 

• FT ' °- 211 > T 5- =0 ' 789 < 53 > 

where r-j<r<r 2 is the tensile region. If r Q defines the crack tip as shown 
in Figure 4, for r Q >r^ in the inner edge crack and for r Q <r 2 in the outer 
edge crack case the crack tip will be in the tensile zone of the cylinder 
and k will be positive. However, in this case the crack surfaces will 
still be partially closed and, hence, the problem is a crack-contact prob- 
lem with the depth e of the contact zone being an unknown constant. The 
physical condition which accounts for this unknown is the smooth closure 
condition of the crack surfaces at r=a+e or r=b-e. Thus, the problem may 
now be treated as an embedded crack problem with the crack surface trac- 
tion p as given by (50) and c=a+e, d=r Q , (r Q >r-]) for inner edge crack, and 
c=r Q , ( r 0 <r 2 )» d=b-e for outer edge crack case. In each case, the smooth 
closure condition to account for the unknown constant e may then be 
expressed as 

k(a+e) = 0 , k(b-e) = 0 . (54) 

For two cylinders with thickness ratios a/b = 0.7 and 0,9, Figure 5 
shows the stress intensity factor k(r Q ). Note that k 1s pos1tive for 
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r^<r Q <b in the inner edge crack case and for a<r 0 <r 2 in the outer edge 
crack case where r^ and r 2 are given by (53), Also note that initially as 
the crack length £ increases (Figure 4} k increases, goes through a maxi- 
mum and tends to zero as the crack traverses the entire cylinder wall. 
Similar results have been observed for flat plates [18], 

For a/b = 0.9 the length e of the contact region is shown in 
Figure 6. For example, in the case of the internal edge crack e=£, 
that is the crack is fully closed for £<r-j-a where (r^-a)/(b-a) = 0.211. 
For £>r-j-a, e decreases monotonously and tends to zero (asymptotically) 
as £ -*> b-a. 
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Table 1 . Stress intensity factors for a symmetric embedded 
crack in a thick-walled cylinder subjected to axial 
tension and pure bending. ((d-c)/(b-a)=0.1 , a 0 = 
Pco/[iT(b 2 -a 2 )j» a 1 =4Mb/[Tr(b 4 -a 4 )]). 
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0.3 


0.615 

1.020 

0.994 

BBS 


0.662 

0.4 

0.06 

0.67 

1.017 

0.996 



0.711 

o 

• 

U1 

0.05 

0.725 

1.014 

0.998 

BBS 


0.761 

0.6 

0.04 

0.78 

1.012 

1.000 

0.799 

I ' 


0.7 

0.03 

0.835 

1.010 

1.002 

0.851 

■ 


CO 

• 

o 

0.02 

0.89. 

1.009 

1.003 

c 

1.903 

fl 


0.9 

0.01 

0.945 

1.007 

1.005 

c 

1.954 

B 


+1.0 

-*0 

+1 

+1.006 

+1.006 

+1 

.006 

+i 

.006 


Table 2. Stress intensity factors for a symmetric embedded 
crack in a thick-walled cylinder subjected to axial 
tension and pure bending. ((d-c)/(b-a)=0.5, a 0 = 
P»/D r(b2-a 2 )], a 1 =4Mb/[Tr(b 4 -a 4 )]). 


a 

b 

s 

H 

M0B 



HKB30BHH 


EP£3 

cose 


0.1 

0.45 

B 

wm 

1.117 

msm 

mSm 

KB 

0.40 

KB 

in 

1.124 

■H 

— 

H 

0.35 

0.475 

1.294 

1.131 


mSM 

m 

0.30 

0.55 

1.268 

1.139 

0.758 


0.5 

0.25 

0.625 

1.247 

1.147 

0.836 

Hi 

0.6 

0.20 

0.70 

1.231 

1.155 

0.911 

0.971 


0.15 

0.775 

1.217 

1.162 

0.984 

1.026 

0.8 

0.10 

0.85 

1.206 



■ 

0.9 

0.05 

0.925 

1.196 

1.178 




+0 

+1 

+1.187 

+1.187 

+1.187 

+1.187 
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Table 3. Stress intensity factors for a synmetric embedded 
crack in a thick-walled cylinder subjected to axial 
tension. ((d-c)/(b-a)=0. 9, a^Pa/OCb^a 4 )]). 


a 

b 

d-c 

“5“ 

c 

F 

r~ i-(- c ) 

"o *?■ 

“»¥ 

0.1 

0,81 

0.145 

3.740 

2.412 

0.2 

0.72 

0.24 

3.300 

2.404 

0.3 

0.63 

0.335 

3.081 

2.410 

0.4 

0.54 

.0.43 

2.945 

2.422 

0.5 

0.45 

0.525 

2.850 

2.437 

0.6 

0.36 

0.62 

| 

2.778 

2.456 

0.7 

0.27 

I 0.715 

2.720 

2.479 

0.8 

0.18 

1 0.81 

2.671 

2.506 

0.9 

0.09 

0.905 

2.627 

2.536 

+1.0 

+0 

+1 

+2.585 

+2.585 


Table 4. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.1, h=b-a» a 0 =P oo /[Tr(b 2 -a 2 )], a^Mb/lXb^-a 4 )]) 


a 

F . 


HHDOSMi 

wmmL Gi 

u 6(0 

Oq/I 

a-j/^ cose 


1 — v ha-jcose 

0 

0.637 


0.128 

0 

0.1 

0.842 

0.123 

0.212 

0.021 

0.2 

0.940 

0.225 

0.244 

0.050 


1.000 

0.334 

0.261 

0.080 

0.4 

1.042 

0.447 

0.272 

o.m 

0.5 

1.073 

0.563 

0.281 

0.142 

0.6 

1.097 

0.680 

0.287 

0.174 

0.7 

1.119 

0.800 

0.292 


0.8 

! 1.138 

0.922 

0.297 


0.9 

1.158 

1.048 

0.302 

0.273 

mam 

+1.189 

+1.189 

+0.310 

+0.310 
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Table 5. Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h«0.1, h=b-a, a 0 =P»/[ir(b 2 -a 2 )], a 1 =4Mb/[ 1 r(b 4 -a 4 )]). 


a 

F 


wmiQmm 


y 5(d) 

KQm 

a-j/jf C0S8 

1-v hr-j cose 

0 

msm 

ID 

■SB 

0.308 

0.1 

BB 

■RIB 

BK 

0.306 

0.2 

■KB 

1.153 

0.302 

0.304 


1.170 

1.149 

0.302 

0.303 

0.4 

1.168 

1.147 

. 0.302 

0.302 

mm 

1.167 

1.147 

0.302 

0.301 

0.6 

1.167 

1.149 

0.303 

0.301 

0.7 

1.168 

1.152 

0.303 

0.302. 

0.8 

1.169 

1.158 

0.304 


0.9 

1.173 

1.166 

0.306 

m BVPV H 

VI. 0 

VI. 189 

+1 . 189 

-*0.310 

-*0.310 


Table 6. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.2, h=b-a, a 0 =P„/l>(b 2 -a 2 )], cr 1 »4Mb/[ir(b 4 -a 4 )]). 


a 

F 



mum 

y 5(c) 



1-v hcr-j cos 

0 

IB 

0.085 


0 

0.1 

■B 

0.153 

B 

0.036 

mm 

■SB 

0.241 

0.441 

0.089 

0.3 

US 

0.342 

0.488 

0.148 

0.4 

1.003 

0.452 

0.524 

0.212 

0.5 

1.055 

0.571 

0.555 

0.280 

0.6 

1.104 

0.699 

0.582 

0.353 

0.7 

1.150 

0.833 

0.608 

0.429 

0.8 

1.198 

0.978 

0.635 

0.511 


1.253 

1.139 

0.666 

0.602 

R0K3I 

->■1.367 

-*-1.367 

-*0.732 

-*0.732 
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• Table 7. Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.2, h=b-a, a^Pc/O^-a 2 }], a-,=4Mb/[>(b 4 -a 4 )j). 


a 

b 

I 

i 

k(c) 

1 

i wmmmom 

V* 

a-j /r cose 

1 -v ha-j cose 

0 

— 

1.314 


0.705 

0.1 

mm 

1.279 


0.685 

0.2 

1.235 

1.253 

1 

0.669 

0.3 

1.231 

1.234 

0.626 

0.658 

0.4 

1.230 

1.222 

0.630 

0.651 

0.5 

1.232 

1.217 

0.635 

0.647 

0.6 

1.238 

1.218 

0.643 

0.648 

0.7 

1.247 

1.226 

0.652 

0.652 

0.8 

1.261 

1.243 

0.663 

0.661 

0.9 

1.285 

1.274 

0.680 

0.677 

+1.0 

+1.367 

+1.367 

+0.732 

+0.732 


Table 8. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.3, h=b-a, a^P^U^-a 2 )], cr-, =4Mb/ [-n-Cb 4 -a 4 ) ] ) . . 


a 

F 

wmmwm 



IHHQHIBOSl 


ff'j/jT cose 

l-v h a Q 

l-v ha-j cose 

0 

BB 

0.127 

I V .; 

0 

0.1 

■131 

0.188 


0.050 

0.2 

0.840 

0.266 

■ 

0.123 

0.3 

0.918 

0.359 

0.707 

0.210 

0.4 

0.991 

0.466 

0.776 

0.308 

0.5 

1.060 

0.587 

0.841 

0.419 

0.6 

1.130 

0.724 

0.906 

0.543 

0.7 

1.203 

0.876 

0.973 

0.681 

0.8 

1.286 

1.053 

1.050 

0.841 

0.9 

1.392 

1.267 

1.150 

. 1.038 

+1.0 

+1.660 

+1.660 

+1.410 

+1.410 
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Table 9 . Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.3, h-b-a, a 0 =P oo /[ 1 r(b 2 -a 2 )], a.,=4Mb/[>(b 4 -a 4 )]). 


a 

b 

HQH 



u sRT 

a o /T 

cr-jyT cose 

1 -v ha-j COS0 

0 

mgm 


0.987 

1.307 

0.1 

■UK 

1.450 

0.976 

1.225 

0.2 

1.328 

1.431 

0.975 

1.165 

0.3 

1.316 

1.381 

0.983 

1.122 

0.4 

1.313 

1.347 

0.996 

1.093 

0.5 

1.317 

1.327 

1.015 

1.078 

0.6 

1.329 

1.323 

1.040 

1.077 

0.7 

1.350 

1.334 

1.073 

1.089 

0.8 

1.384 

1.365 

1.117 

1.120 

0.9 

1.442 

1.428 

1.185 

1.181 

->■1.0 

->■1.660 

-*■1.660 

->■1.410 

-*•1.410 


Table 10. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.4, h=b-a, a 0 =P«,/O(b 2 -a 2 )], a ] =4Mb/[7r(b 4 -a 4 )]) . 


a 

b 


mmomm 


nrasi 

a 0 A 

a-j/^cose 

l-v ha-j cos0 

0 

0.665 

0.171 

0.531 

0 

0.1 

0.754 

0.226 

0.686 

0.062 

ESI 

0.838 

0.296 

0.817 

0.155 

rm 

0.920 

0.383 

0.935 

0.269 

eh 

1.001 

0.487 

1.046 

0.404 

0.5 

1.085 

0.611 

1.158 

0.563 

0.6 

1.174 

0.757 

1.277 

0.752 

0.7 

1.275 

0.928 

1.412 

0.973 

0.8 

1.397 

1.141 

1.580 

1.254 

0.9 

1.568 

1.426 

1.821 

1.638 

Emm 

-*2.112 

-*2.112 

-*2.614 

-*2.614 
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Table 11. Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
(£/h»0.4, h-b-a, cx 0 =Pe»(l>(b 2 -a 2 )], cr^Mb/fr^-a 4 )]) . 


a 

F 

H0M 

1 

1 


HV9HH3G9JBI 

a o& 

a-j/^ cose 

1 -v ha-j COS0 

0 


2.077 

WIHHH 

2.337 

0.1 

1.513 

1.865 


2.078 

EbT 

1.465 

1.715 


1.899 

jm 

1.437 

1,606 

1.392 

1.774 

llli 

1.425 

1.531 

1.419 

1.692 

0.5 

1.427 

1.486 

1.459 

1,646 

0.6 

1.443 

1.467 

1.516 

1.633 

0.7 

1.475 

1.476 

1.594 

1.658 

0.8 

1.533 

1.520 

1.706 

1.730 

0.9 

1.641 

1.626 

1 .891 

1.890 

+1.0 

+2.112 

+2.112 

+2.614 

+2.614 


Table 12. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0. 5, h=b-a, a 0 =P«/|>(b 2 -a 2 )], a-,=4Mb/[-n-(b 1+ -a 4 ) ] ) . 


a 

F 

■OH 


— 

1 

1 


a- j/j cos 0 

1 -v ha-j COS0 

0 

0.691' 

m ■ 


0 

0.1 

0.775 



0.074 

0.2 

0.859 



0.187 


0.945 

■Ear, 


0.329 

0.4 

1.035 

0.517 


0.503 

0.5 

1.131 

0.643 

1.523 

0.718 

0.6 

1.239 

0.799 

1.718 

0.985 

0.7 

1.366 

0.989 

1.954 

1.318 

0.8 

1.529 

1.243 


1.775 

0.9 

1.779 



2.472 

+1.0 

+2.826 


+4.950 

+4.950 
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Table 13. Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.5, h=b-a, a^OC^-a 2 )], ^Mb/CTt^-a 4 )]) 


a 

F 

1Q1 

mmamm 

~ I ' 

mamms-m 

■a 

cr-j/T C0S8 

1 — v ha-j cose 

0 

mm 

2.929 

1.989 

4.291 

0.1 

mm 

2.468 

1.901 

3.548 

0.2 

1.667 

2.163 

1.873 

3.076 

0.3 

1.610 

1.946 

1.884 

2.753 

0.4 

1.580 

1.780 

1.924 

2.546 

0.5 

1.572 

1.707 

1.994 

2.426 

0.6 

1.586 

1.660 

2.097 

2.382 

0.7 

msm 

1.658 

2.248 

2.417 

0.8 

BEB 

1.712 

2.478 

2.558 

0.9 

1.878 

1.867 

2.888 • 

2.904 

-►1.0 

->2.826 

->2.826 

->4.950 

->4,950 


Table 14. Stress intensity factors and crack opening displace- 
ments for an internal edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending. 
U/h=0.6, h=b-a, a 0 =P»/[ir(b 2 -a 2 )], a-,=4Mb/[>(b 4 -a 4 )]) 


a 

F 

■IQB 

■ncim 


1 

1 

a 0 /r 

cr-j/s, cose 

1-v ha-j cose 

0 

mm 

0.265 

0.877 

0 

0.1 

H 

0.314 

1.077 

0.087 

0.2 

m 

0.378 

T. 283 

0.222 

0.3 

1.000 

0.459 

1.493 

0.393 

0.4 

1.099 

0.561 

1.716 

0.611 

0.5 

1.208 

0.691 

1.966 

0.890 

0.6 

1.333 

0.856 

2.261 

1.253 

0.7 

1.484 

■ 

2.636 

1 .729 

0.8 

1.688 

1 1 E . 9 1 

3.166 

2.431 

0.9 

2.025 

1.824 

4.091 

3.623 

EH9 

->4.035 

->4.035 

->9.965 

-9.965 


- 23 - 
























Table 15. 


Stress intensity factors and crack opening displace- 
ments for an external edge crack in a thick-walled 
cylinder subjected to axial tension and pure bending, 
U/h=0.6, h=b-a, a^P-ZO^-a 2 )], a 1 =4Mb/I>(b‘*-a 4 )]) 


a 

F 


HOHi 


y 6(d) 


U-j/i" COS0 

1-v ha-| cos9 

0 

2.478 

4.579 

2.798 

8.542 

0.1 

2.159 

3.527 

2.589 

6.365 


1.977 

2.880 

2.510 

5.095 

0.3 

1.866 

2.460 

2.506 

4.307 

0.4 

1.802 

2.194 

2.557 

3.833 

0.5 

1.773 

2.021 

2.661 

3.550 

0.6 

1.776 

1.923 

2.826 

3.424 

0.7 

1.818 

1.895 

3.080 

3.455 

0.8 

[ 1.918 

1.949 

3.491 

3.692 

0.9 

2.153 

2.155 

4.289 

4.354 

VI. 0 

-4.035 

-4.035 

->9 . 965 

->9.965 
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Table 16. The effect of Poisson's ratio on stress intensity 
factors and 6(c) when loading is uniform tension 
(a/b=0.5, */h=0.3, 0 o =P«/[ir(b 2 -a 2 )], h=b-a). 



v=0 

v=0. 1 

v=0.2 

CO 

• 

o 

II 

7 > 

v=0.4 

v=0.5 


1.048 

1.051 

1.055 

1.060 

1.067 

1.076 


0.814 

0.821 

0.831 

0.841 

0.854 

0.870 


Table 17. The effect of Poisson's ratio on stress intensity 
factors and 6(c) when loading is pure bending 
(a/b=0.5, z/h=0.3, o^Mb/CirM-a 4 )], h=b-a). 



v=0 

v=0. 1 

EZI 

v=0.3 

v=0.4 

v=0.5 

k «n 

0-j /£ COS0 

0.574 

0.577 

0.582 

0.587 

0.594 

0.602 

-JL... 5 (P) 
1-v hcr-jCose 

0.395 

0.402 


0.419 

0.430 

0.443 


- 25 - 































APPENDIX A 

The functions_C in (s) which appear in (28) (m^Cs)) = M(s) = F TCs) » 
F(s) = (f..(s)), (i,j) = 1 ,6. The coefficients f.j of the matrix F are 

fj-j = rs I n (sr) - (2-2v+n) I n+1 (sr) (A-l) 

f i2 = -[rs ysr) + (2-2v+n) K^sr)] (A-2) 

f i3 * s w sr) + f y sr) (A ' 3) 

f i4’ r K n (sr) ‘ s K n+l (sr) (A ‘ 4) 

f 15 - - 7 (A ' 5) 

f 16 = ' r l V sr > ■ 1 = !’ 2 > < A " 6 ) 

for i = 1, r = a and for i = 2, r = b. 

f i i = - ( - n y - ) (4-4v+n) I n+1 (sr) + (2-2v+n)s I n (sr) (A-7) 

f i2 = " ■^T^' (4-4v+n)K n+1 (sr) - (2-2v+n)s K n (sr) (A-8) 

(n ' ,) V 51 -’ (A ‘ 9) 

f 14 * -f K n+ -j( sr ) + fr (n-1) K„(sr) (A-10) 

f 15 - - Q (n-l)+s 2 )I„(sp) + 25. I„ +1 (sr) (A-ll) 

f 16 - - (fg- (n-l)+s2)K n (sr) - 2l K n+ ,(sr) (A-12) 
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for i = 3,4; when i = 3, r = a and when i = 4, r = b. 

f n = (3-2v)s I n (sr) - (4-4v+n) + rs 2 ]I n+1 (sr) 

(A-13) 

f i2 = - {(3-2v)s ^(sr) + [-£U±D- (4-4v+n) + rs 2 ]!^^)} 



(A-14) 

f t3 = - tC^+s^tsr) -f I n+ ,(sr)} 

(A— 15) 

f 14 • - U 2 ^ 0 - + s 2 ]K„(sr) + f K n+ ,(sr)} 

(A-16) 

f 15 = f C^I„(sr) +S I„ +1 (sr)] 

(A-17) 

fl6 . |l [Indl ysr) . s ^drj] 

(A— 18) 


for i = 5,6; when i = 5, r = a and when i =6, r = b. 

The terms Gj (s, t) in (28) are obtained by evaluating the Bessel 
integrals which result from the substitution of (27) into equa- 
tions (8), (11), (12), (14) and (21 )- (24) . These definite Bessel 
integrals are evaluated by differentiating a related integral 
given in [17], that is, by using 

| a v-„+2n+l . (-1) V-" +2n yas)K v (ts) 

0 

a > 0, Re s > 0, a < t < «, 

Re y - 2n+l > Re v > -n-1 , n = integer. (A-19) 
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Taking the derivative with respect to s on both sides, the 
integral becomes, 

ct> p 

j a V ^ +2n+ ^ (s 2 +a 2 ) 2 J^( act) (ta)da = — {(v-p+2n) * 

o 

* s v_p+2n “ 2 ^ (as ) Ky ( ts ) + s 'w-i i +2n-T[a ^ +l ( sa) + 1L r(, sa )> 

* K y (ts) + s v “ y+2n ‘ 1 I y (as)[j K v (st) - t K v+1 (st)]} 

a > 0 ,Re s > 0, a < t < «, (A-20) 

Re y-2n+4 >Re v > -n-1 

Thus, the terms Gj(s,t) (j=l,6) are expressed as: 

G-j (s,t) = K n (st){ 2 ^ ^ - I n (as) + (4+3n)s I n+1 (as) 

+ s 2 a I n+ 2^ as )> " t K n +l^ st) {I eF ! n (asJ + s2 Vl^* 

(A— 21 ) 

G 2 (s,t) = I n (st){^In±IlK n (bs) - (4+3n)s K n+1 (bs) 

+ s 2 b K n+2 (bs)} + t I n+1 (st) K n (bs) - s 2 K n+1 (bs)} 

(A-22) 

G 3 (s,t) = £ {K n (st)[^lllvj^ yas) + s (2v+3n+l )I n+1 (as) 

+ s 2 a I n+2 (as)] - t Kn +1 (st)[^=^- I n (as) + s 2 I n+1 (as)]> 

(A-23) 
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S 4 (s,t) =■ £ tI n (st) C^g^ -K n (bs) - s(2v+3n+l)K n+1 (bs) 

+ s 2 b K n+2 (bs)] + t K n (bs) - s 2 K n+1 (bs)]} 

(A-24) 

S 5 (s,t) = K n (st){-rMBjTKni2l + 2(n+1)s 2 ]I n (sa) 

+ if- (-n 2 +3n+2(v+l))-s 3 a]I n+1 (sa) + s 2 I n+2 (sa)} 

-t I n (as) + ^ I„ +1 (as)} (A-25) 

G 6 (s,t) = I n (st){- [ 2n ( n -^ n ^) + 2(n+1)s 2 ]K n (bs) 

+ [f (n 2 -3n-2(v+l)) + s^b^tbs) + s 2 K n+2 (bs)} 

+ t I n+ ,(st){-[2^- + s 2 ]s ybs) - £ K n+1 (bs)} ( A— 26 ) 




APPENDIX B 


The functions v-(s,t) which appear in (30) 

J 

v j = Y jK Q K for j = 1 » 3 ’ 5 * K = 1,2 (B-1). 

Vj = Yjk Q k+2 for j * 2,4,6 , K=l,2 . (B-2) 

For example, 

V 1 = Y Vl q l + Y 12 Q 2 
v 2 = t 21 Q 3 + y 22 q 4 
v 3 = Y 31 Q 1 + y 32 Q 2 • etc ' 

Y 11 = x n (as ) + (4+3n)s I n+1 (as) + s 2 a I n+2 (as) (B-3) 

>12 a fV B l ts!l riW (*- 4 > 

y 21 = 2n b n+1 ? K n (bs) ' (4+3n)s l W bs) + s2b K n+2 (bs) (B ' 5) 
Y 22 =! T K n (bs) • s2K n+1 (bs > (B-6) 

y 31 = ? I n (as) + s(2v+3n+l)I n+ ,(as) + s 2 a I„ +2 (as)] 

(B-7) 

y 32 * a I„(as) + s 2 I n+1 (as)] (B-8) 

y 41 * F [ 2(n 'b )(y+n) K n ( bs > - s(2v+3n+l )K„ +1 (bs) + s 2 b ^(bs)] 

(B-9) 
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I42 * F V bs) - * 2 W bs >3 < B - 10 ) 

y 51 - .[ ialaiTKlttvl ♦ 2(n+l)s 2 ]I n (sa) + [f (-n 2 +3n+2(v+l)) 

- s 3 a]I n+1 (sa) + s 2 I n+2 (sa) (B-ll) 

y 52 - -C 2 ^ + S*]s I n (as) ♦ f I n+1 (as) (M2) 

y 61 * + 2(n+l)s 2 ]K n (bs) + [f (n 2 -3n-2(v+l)) 

+ 5 3 b]K n+1 (bs) + s 2 K n+2 (bsj (B— 13) 

y 6 2 * -c 2 ^ + s 2 ]s K^bs). - £ K„ +1 (bs) CB-14) 

The terms in (B-l) and (B-2) denoted by are defined as, 

Q-|(s,t) = - i [nt - It K n _ 1 (st)S 1 >n (1s-t)3 

(B-15) 

Qj(s.t) * JT [(n+2)t K n+1 (st)S l n (1st) - It K n (st)S 2 n+1 (1st)] 

(8-16) 

Q 3 (s,t) = - j [nt I„{st)S 0>n _ 1 (1st) + it I n _ 1 (st)S 1>n (ist)] 

(8-17) 

Q 4 Cs,t) » - Jr [(n+2)t I n+1 (st)S l n (ist) + It !„ (st)S 2 n+1 ( 1st)3 . 

(B-18) 

The Lomrnel functions of imaginary argument which are encountered 
in Q k are not in a form which is convenient for numerical 
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computation, but with the aid of certain identities involving 
the Struve function Q K may be expressed in 


terms of real functions with real arguments. The terms Q K 


K=l,...,4 for five harmonics n=0,...,4 are given 

below. 

n=Q 


Q-j (s ,t) = K-j(st) 

(B-19) 

Q 2 (s,t) -fl^Cst) +~K 0 (st) 

(B-20) 

Q 3 (s,t) = j J i (st) 

(B-21 ) 

Q 4 (s,t) = - f| I-j Cst) + ^ I 0 (st) 

(B-22) 

1 

let: u Q (z) =z | e~ zp (l-p 2 )"^ 2 dp 

0 

(B-23) 

U-, (z) = | e“ zp (1-p 2 ) 1/2 dp 
0 

(B-24) 

u 2 (z) = z 2 | e" zp (1-p 2 )^ 2 dp, 
0 

(B-25) 

then for n=l : 


(s,t) = ^ [K-j (st)u Q (st) + K Q (st)u^ (st) J 

(B-26) 

Q 2 ( s , t ) = [Si^CstJu-j (st) + K.j(st)u 2 (st)] 

(B-27) 
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(B-28) 


Q 3 (.Sjt) = I Cl^stju^st) - I 1 (st)u 0 (st)] 

Q 4 (s,t) = p- [SlgCstJu^st) - I 1 (st)u 2 (st)] (B-29) 

n=2 

Q-j(s,t) = Kgtst) - j D - -[ilF ]K l (st) (B ’ 30) 

Q 2 (s, t) = gT Cl " XstF^ K 3( st ) + s 7 ^- st " st + (st^’V 5 ^ 

(B-31) 

Q 3 (s , t) = - Jr I 2 (st) + | Cl - j 5 |yzGi 1 (St) (B— 32) 

Q 4 (s»t) = ^ [1 - -[i£p-]I 3 (st) ■+ J [st - J + 

(B-33) 

n=3 

Q-j (s,t) = I' {K 3 (st)[3u 0 (st) - J u-j (st)] + K 2 (st)[3u 1 (st)- u 2 (st)]} 

(B-34) 

Q 2 (s.t) = jr {K 4 (st)[|^ u 2 (st) - 15 u 1 (st)3 

+ K 3 (st)[st-15 u Q (st) + -^r u-j (st) - -^|p- u 2 (st)]} (B-35) 

Q 3 (s,t) = | {I 3 (st)[3 u Q (st) - J ^ (st)] 

- I 2 (st)[3u-,(st) - ~ u 2 (st)]} (B-36) 

Q 4 (s»t) = {I 4 (st)[|| u 2 (st) - ISu^st)] 

+ I 3 (st)[-st+15u 0 (st) - ™ u-,-(st)- + u 2 (st)]} (B-37) 
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n=4 

Tt + Ti§F ] M s t) - ! [1 - THF 

+ Xi?T r - lK 3 (st) (B " 38) 

Q 2 (s,t) “ f| [1 - ^fpr + ^If^Kjtst) + pr [st - || 

+ tSf ■ fnF ]K 4 (st) < B ' 39 ) 

Q 3 (s,t) .f -Jp, 

+ iilf ,r]I 3 (st) < B - 4 °) 

Q 4 ( S .t) =^[1 - Ti if r + T ||f T ]I 5 (st) 
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APPENDIX C 

Evaluation of LommeVs Functions 

In the problem considered, LonmeVs functions S y#v arise from the 
indefinite integrals of the form 

| z w J v (z)dz = (u+v-1)z J v (z) s M -l,v-l (2 ) - z J v-l (z)s v,v (z) 

(C-l) 

t 

S^z) = z u-l n ,lBd^i t [(v-l) z -^[(v-3) z -v 2 ] . ...j 

(C-2) 


Referring to [12] it can be shown that the Lornmel's functions of real 
argument may be expressed as 


n = 0 


n = 1 



1,0 (z) = 1 

(C-3) 

0>0 (z) - \ (i4t*r * zt dt 

(C-4) 

i, i(z) = z f (1+t 2 ) 2 e“ zt dt 

(C-5 ) 

* ^ 0 


o 

N 

II 

IM 

1 

(C-6 ) 

I 1>2 (z) = 1 + 4z’ 2 

(C-7 ) 
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f S o,2^ = 2z 1 S 1,l( z ) " S o,o^ z ^ ^ C ' 8 ^ 

n= 3 

1 s 1>3 ( z ) = 4 + (24 z- 2 -.3)S 1#1 (z)-12z- 1 S oo Cz) 

(C-9 ) 

r s o,3^ z ^ = Z_1 + 8z " 3 (c-io) 

n=4 < 

( S lj4 (z) = 1 + 16z“ 2 + 192z“ 4 (C-ll) 

i 

S 0j4 (z) = 8 z“ 1 +(48z" 3 -8z' 1 )S 1j1 (z)+(1-24z‘ 2 )S 0>0 (z) 

n=5 i ( C- 12) 

1 S 1 ,5^ z ) * -4+320 z _2 +(1920z‘ 4 -360z‘ 2 +5)S 1 ^z) 

+ (60z" 1 -960z‘ 3 )S OjO (z) (C-13) 


The details for evaluating S asymptotically for small and large argu- 

u > v 

ments may be found in [19]. 
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APPENDIX D 

Evaluation of the Integral in (32) 

As noted in Section 2, to evaluate the integral in (32) it is neces- 
sary to analyze the integrand asymptotically. By using the asymptotic 
expansions for J n (z), S (z), S( x )> and c ( x ) as 9 iven by [12-14] and by 
evaluating the integral in A<a<» asymptotically, (32) may be expressed as 


0 

- J n _ 1 (at)S 1 n (at)]J n (ar)- — ! — [C(i^t)(sin(ar)+(-l) n cos(ar)) 
* /2irra 


- S(v^Tt)(cos(ar)+(-l) n sin(ar))]}da + [A-j+Ag+Ag], n=0,l,... 

(D-l) 


where A is a "large" number and A-j , Ag, and Ag are given by 


n /t 


A, ■ 7 (?) {-CKA(t-r)) + (-l) n [f - Si(A(t+r))] 

f 4 n£-l J_ _ ( 4 n 2 - 1 ) ( 4 n 2 - 9 ) ,1 + 1 ,, * 

L 8 tr T28 + friJ 

* [ C ps^t-r) + (-T) n A(t-r) 

+ (-1) n ^I cos^t+ri + c , (A(t . r)) 

. t-D" [| . si(A(t+r))]] + (1 - 1)* 

* C Al J- ACt-r) , + ( . 1} n cos _A (t*r) . (t . r)c1(A(t . r)) 
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- (-l)"(t+r) [f - SI (A(t+r))]]> (D-2) 



= 1 (h H <r- (4n 2 -1)(4n 2 -9) -2 

7T V 1L 128 r 


(4(n-1) 2 -l)(4(n-l ) 2 -9)-64n 2 t -2-, r sin A(t-r) 
128 A 


- (t-r)Ci(A(t-r)) - (-1)" cos^tjr]. + ( .,)n (t+r)( j 

- Si(A(t+r)))] + 4(n-1) 2 -l t -l [_ci' (A(«-r) ) 


+ (-1)" (f- Si(A(t+r)))] + ^lr' 1 [Ci (A(t-r)) 

+ (-l) n (f - S1(A(t+r)))] + (4n 2 -1?(4|n-l) 2 -1) (rt) -1, 

*[ s1n * (t ~r) . ( t -r)ci(A(t-r)) + ( -i) n c - ° - s gted 

-(-l) n (t+r)(f - S1(A(t+r)))]} (D-3) 

A 3 = - [C(/Sr5 - S{/SF)3 - i (rtf’* * 

- S1(A(t+r))] - CiCA(t-r))] - i r ' H t ' % * 

*[- sin A(t-r) + (t . r)c1(A(t . r)) + H) n cos A(t+r) 

- (-1)" (t+r)(|- S1(A(t+r)))] , n = 0,1,2 CD-4) 


The definite integral in (D-l) is evaluated by using Gauss-Legendre 
quadrature with the upper limit A = 200. 
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z 



Fig. 2 Stress intensity factors for concentric ring 
shaped cracks in an infinite medium subjected 
to axial extension (n=0) and pure bending (n=l). 
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Fig. 4 Geometry of edge cracks in a cylinder wall sub- 
jected to residual stresses. 



k(r 0 ) 



b-a 

Fig. 5 Stress intensity factors for edge cracks in a hollow 
cylinder subjected to residual stresses, a/b = 0.9, 
a/b = 0.7. 
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b-a 


Fig. 6 Crack contact length e for a hollow cylinder under 
residual stresses with either internal or external 
circumferential edge crack, a/b =0.9. 
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